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Abstract. We apply a 3-dimensional approach to describe a new parametrization of 
the L-operators for the 2-dimensional Bazhanov-Stroganov (BS) integrable spin model 
related to the chiral Potts model. This parametrization is based on the solution of the 
associated classical discrete integrable system. Using a 3-dimensional vertex satisfying 
a modified tetrahedron equation, we construct an operator which generalizes the BS 
quantum intertwining matrix S. This operator describes the isospectral deformations 
of the integrable BS model. 

PACS numbers: 05.45-a, 05.50+q 
Introduction 

The aim of this paper is to describe in detail the interrelation between specific 3- 
dimensional (3D) and 2-dimensional (2D) integrable lattice spin models. In particular, 
we shall make use of 3D techniques in order to derive the isospectral transformations 
of a 2D model, which would be difficult to find directly. The 3D model we shall exploit 
is the generalized Zamolodchikov-Bazhanov-Baxter (ZBB) model PI El 0] in the vertex 
formulation [U EJ E] ■ This model describes chirally interacting Z^v-spins placed on 
the links of a 3D cubic lattice. The corresponding 2D model will be the integrable 
Bazhanov-Stroganov (BS) model [H], which is related to the integrable Chiral Potts 
(CP) model HHHUHg. 

We shall use the approach to the generalized ZBB model developed in E] . The 
dynamical variables are affine Weyl pairs u, , Wj which live on the links j of an oriented 
3D lattice. The cornerstone of the approach is the explicit construction of a canonical 
map of the triplet of Weyl pairs on the three incoming links of a vertex to the Weyl 
pair triplet on the three outgoing links. This map defines the Boltzmann weights and 
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by construction satisfies the tetrahedron equation. The canonical map is uniquely 
determined postulating a Baxter Z-invariance and a specific linear problem for the 
Weyl variables. 

The generalized ZBB model with Z^-spins is obtained taking the Weyl variables 
to commute to the iV-th root of unity. Then the Weyl centers and are classical 
dynamical variables of a classical discrete integrable system of Hirota form which is 
determined by the canonical mapping and boundary conditions. This system can be 
solved using standard tools of algebraic geometry, i.e. G functions and Fay-identities. 
Using the rational limit of the B-functions this is handled in a practical explicit way. 

The Weyl operators at iV-th root of unity are represented by N x N matrices. 
Then the canonical mapping decomposes into a functional mapping of the centers and a 
finite-dimensional transformation given by iV 3 x iV 3 matrix. It was shown in jH] that the 
matrix element of this matrix coincide with the Boltzmann weights of ZBB model. In 
this approach the ZBB model appears in the special case that the functional mapping 
is trivial, i.e. that trivial solutions to the Hirota-type equations are chosen. Non-trivial 
solutions for the functional mapping mean non-trivial classical dynamics of the Weyl 
centers, in particular, solitonic solutions. The final aim is to find separated variables. 

It is well-known that the few-layer ZBB model, with quasiperiodic boundary 
conditions in the direction orthogonal to the layers, is related to the integrable chiral 
Potts model. The main aim of this paper is the explicit construction of various properties 
of the Bazhanov-Stroganov quantum chain directly from the linear problem and the 
canonical mapping operator of the generalized ZBB model. The linear problem leads 
to the BS L-operator. The quantum intertwining operator of the BS model is obtained 
as the product of two 3D canonical mapping operators. In case of the trivial functional 
mapping this gives the well-known BS S-matrix. However, this is generalized if non- 
trivial classical dynamics is taken into account. Intertwining through the whole BS chain 
leads to isospectrality transforms of the transfer matrix. A special case of the BS model 
is the relativistic Toda chain, for which isospectral transforms have been constructed 
already in ^H]- An important advantage of the 3D approach to 2D problems is the 
flexibility regarding the choice of parametrization. The CP parametrization turns out 
to be less convenient for the dynamical case than a parametrization using simple cross- 
ratios and rational B-functions. 

The paper is organized as follows: In section^ we summarize the main features and 
formulae of the models considered. Then in section 121 the L-operator and the quantum 
intertwining relation for the BS-model will be derived using the canonical mapping 
approach to the ZBB-model. A new parametrization of the BS-intertwining matrix in 
terms of cross-ratios is introduced. In the following section |3] we introduce a classical 
counterpart of the BS-model and find the transformation realizing the intertwining of 
two Lax-operators, using the functional mapping of the 3D vertex ZBB-model. Section 
0] starts stating the main result of the paper, the explicit formula for the isospectral 
transformation of the BS-model. The proof of this proposition is given in the following 
subsections. Section El summarizes the results. 
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We start with a summary of some basic features of the models considered in the later 
sections. This will also serve to establish the notation. 

1.1. Vertex formulation of the generalized ZBB-model 

In the vertex formulation of the ZBB-model jl] the quantum variables are attached to 
the links j of a 3D oriented lattice. They are taken to be elements (u.j,Wj) of an 
ultra-local affine Weyl algebra at root of unity: 



Uj- • Wj = u Wj ■ Uj] u N = l; NeZ; N>2 (1) 

and Uj • Wj = Wj • Uj for % ^ j. We also attach a scalar variable Kj to each link j and 
denote these variables together as 




Figure 1. Left: The six links of the basic lattice intersecting in the vertex A, 
intersected by auxiliary planes (shaded) in two different positions: first passing through 
roi, rt>2, r03 and second through tr^, tt> 2 , tt> 3 . The second position is obtained from the 
first by moving the auxiliary plane parallel through the vertex A. The Weyl variables, 
elements of n>i, ro-, live on the links of the basic lattice. 72-123 can be considered to 
be attached to the vertex A, it maps the left auxiliary triangle onto the upper right 
one. Right: the auxiliary plane in the neighborhood of tt>i, showing the "co-currents" 
($ a |, ... , ($(2 1 in the four sectors cut out by the directed lines W2W1 and W3W1. 
The Linear Problem relates these four adjacent co-currents according to the values of 

Wi = (ui,Wi,Kl). 

Fig. ^ shows on the left the three Weyl variables to j on the ingoing links of a vertex A 
and the three variables ro^ on the corresponding outgoing links. 

In the approach of [S] the basic object of the generalized vertex ZBB-model is the 
operator 7^123 mapping canonically the triple affine Weyl algebra on the ingoing links 
to the corresponding triple Weyl algebra on the outgoing links. This mapping is an 
invertible rational mapping: For any rational function \1/ of the ui, . . . , W3, we define 

(7^123 *) (ui,Wi,U 2 , ••• ,W 3 ) = #(ui, wi,U 2 , . . ., W3). (3) 
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In jSJ |B] , 7^123 has been determined uniquely postulating a Baxter Z-invariance (lines 
may be shifted with respect to each other) and a linear relation ("Linear Problem") 
between the four "co-currents" attached to the four adjacent sectors around each tOj 
in the auxiliary plane. The right hand part of Fig. ^ shows the auxiliary plane in the 
neighborhood of tDi and the corresponding four co-currents. The "Linear Problem" at 
roi is taken to be 

= ($ | +w 1/2 ($ fe | Ui +($ c |Wi +Ki($rf|UiWi, (4) 

analogously at all links j. The lines in the auxiliary plane have a direction (we shall not 
go into the rule here) so that e.g. the co-current appearing in (J3J) multiplied by 
is the one between the arrows. 

For the 2D auxiliary plane the relation (@J) contains analog information as does in 
the standard ID quantum chain case the QIS L-operator relation 

Before giving the explicit formula for the mapping operator, we use the fact that 
for uj a iV-th root of unity, the affine Weyl operators can be represented by N x iV 
matrices. Omitting for a moment the index j (because of the ultralocality the full 
representation space is just a direct product), we write 

u=uX; w=wZ; u, w e C; XZ = wZX, (5) 

and shall use the natural basis 

X|/3) = Z|/3) = \(3 + l); ( a \(3}=5 a ^. (6) 

Clearly = Z N = 1 . The iV-th powers of the Weyl elements are centers and we 
shall denote them by Uj, Wj : 

uf = < ee Uj- wf = wf ee Wj. (7) 

Now (uj+Wj ) N = Uj + Wj , and the mapping TZu3 implies a purely functional mapping 
( f) 

TZi 23 °f the centers Uj, Wj , or taking iV-th roots, of the Uj, Wj : 

(n { / 2 ] 3 o ip\ ( Ul , Wt, u 2 , w 3 ) d = i)(u[, w' x , u' 2 , ... ,w' 3 ). (8) 

The remarkable feature (observed in [23 J arises that 7^123 decomposes into a matrix 
conjugation R 12 3 (this is a iV 3 x A^ 3 -matix) and a purely functional mapping TZ^/ 2 \ '■ 

n 12 3 o * = r 123 (ng\ o r~ 2 \ . (9) 

The matrix R123 can be given compactly in terms of the Bazhanov-Baxter cyclic 
functions w p (n) (not to be confused with w in (jSJ)) which depend on the Z^-variable n 
and on a point p = (x, y) restricted to a Fermat curve: 



w p {n) y v v 



-I 2^+^ = 1; neZ N ; n > 1; w p (0) = 1 . (10) 
w p {n — 1) 1 — uo n x 

The cyclic property w p (n + N) = w p (n) is guaranteed by the Fermat curve restriction 
in (|1U|) . The functions w p (n) are root of unity analogs of g-gamma functions and can 
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be used to develop the theory of the corresponding g-hypergeometric functions (see e.g. 
jUHH]). One can show 000123 that R 12 3 can be written as a weighted cross-ratio 
of four of these cyclic functions. In components: 

T>hhh _ X (n - ll)h W pi {l 2 ~ tl)w p2 (j 2 - 

414213 " d ^ + ^ w P3 (h - h)w P4 (i 2 - j x ) ■ W 

Here Pi = (xi,yi), i = 1,2,3,4 are four points on the Fermat curve (jlOj) which are 
related by the constraint jlj 

Xi X 2 — UJ X 3 X4 . (12) 

These Fermat points can be expressed in terms of the parameters Uj, Wj, k,j, u'j, w'j of 
the initial and final Weyl pairs tOj, xv'j (j = 1, 2, 3) : 

-1/2 U 2 -1/2 «2«2 -1 U 2 

Xi = uj 1 ; x 2 = uj 1 — — ; x 3 = uj — ; 

K\U\ U\ U\ 

y 3 _ KgWx _ y^ _ 1/2 k 3 w 3 

u[ , u' 2 , u 3 N are related to the initial variables by the functional transformation 
Defining Kj = , the mapping TZ^/ 23 is explicitly given by 

U{ Wi _ K 2 U 2 W 2 



, . 1 (13) 

2/i ^3 Vi w 2 



U x W 3 K 1 U 1 W 2 + K 3 U 2 W 3 + K X K 3 XJ X W 3 

W l _W[ W1W3 

Wi~W 2 ~ W X W 2 + U 3 W 2 + K 3 U 3 W 3 ' 
U' 2 U 3 U,U 3 



(14) 



U 2 ~ U! 3 U 2 U 3 + U 2 W X + K X U X W X ' 

If we need the u'j rather than the UL which is the case when we calculate the Fermat 
points via ()13|) . we have to take A^-th roots. The choice of phases is restricted by the 
fact that the complete mapping TZi 23 leaves the following three products invariant [25]: 

W[ w 2 = wj Wj; u 2 u 3 = u' 2 U3; U! w 3 1 = u'jWj 1 . (15) 

Considering a 3D model of Zjv spins on the links of the lattice, the R^H^ 3 can be taken 
to be the (generally not positive) Boltzmann weights of the vertices. Via the Fermat 
parameters, these depend on the scalar parameters u\, u 2 , . . . ,K 3 . Each solution of the 
functional equations gives rise to an integrable 3D model. 

It can be seen, that by construction, the mapping 7Zi 23 satisfies the tetrahedron 
equation and that R123 solves the modified tetrahedron equation, see e.g. [26J. 



1.2. Integrable Chiral Potts model 

The integrable chiral Potts model (CP) is defined on a 2D lattice with Z^r spins Oj 
attached to the vertices. There are two sets of directed rapidity lines p, p', and q, q' 
which cross on the links of the lattice, see Fig. El If the edge linking the spins is between 
the rapidity directions, the pair interaction Boltzmann weight (which depends on both 
rapidities crossing on the link) is W pq (o-j — cry). If the edge is to the right of the rapidity 
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directions, it is W pq (o~j — a^). Such Boltzmann weights which satisfy the star-triangle 
relation (for the explicit form of R pqr see e.g. [TBI ) 

^W qr (d-b)W pr (d-a)W pq (c-d) = R pqr W pq {b- a)W pr {c-b)W qr {c- a) (16) 



P 



V 



T q {a,a f ) 















\W P > q 









w, 



/<</ 



, 1 



Q 



<j\ a 2 cr 3 

Figure 2. The square diagonal directed lattice with the Z^-variables at the vertices 
and the Boltzmann weights W pq on the right pointing links and W pq on the left pointing 
links. The weight corresponding to the line from 02 to a' 2 is W vq {o2 — <j' 2 ), analogous 
for the link from a 2 to a[ it is W v i q {a 2 —&[)■ Dashed lines are the rapidity lines which 
indicate the parameter dependence of the Boltzmann weights. For simplicity, we show 
only the special case of alternating rapidities p,p' and q,q'. There are two different 
transfer matrices T and T. 



have been constructed in ^0] arid are defined by the difference relations (n G Zjv) 

Wpgjn) _ f lh\ Vi ~ uTIx p . W pq {n) = ujx p - u n x q 

W pq (n-1) \n q ) y p -uj n x q W pq {n-1) P q y q - u n y p 

We shall use the normalization W pq (0) = W pq (0) = 1. The parameters appearing in 
(fTTjl are constrained to the high-genus curve 

< + yr = *«y? r + 1 ); kx N q =i-k'^ q N ; k y ? = i-k>^, (is) 

(same for x p , y p , fi p ) where k and k' are fixed temperature-like parameters related by 
k 2 + k' 2 = 1 . The constraints (|T%|) guarantee the cyclic property 

W pq (n +N) = W pq (n); W pq (n + N) = W pq (n) . (19) 

The star-triangle relations |T6|) are quite special since the three rapidities involved 
appear not as differences as usual, but each separately. Due to this feature, many 
standard techniques cannot be applied straightforwardly to the CP model. Functional 
relations involving the BS-model discussed in the following have been crucial for 
obtaining analytic solutions for the CP-model, see e.g. 0111112011211122] and references 
therein. 
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1.3. Bazhanov-Stroganov model 



R{\v) 



(20) 



The CP model has been found to be intimately related to the six-vertex model in a 
seminal paper by Bazhanov and Stroganov [£.:. They first noticed that the twisted 
six-vertex R-matrix 

( X-uv 0\ 

u)(\-v) A(l-w) 

-u) X-v 

\ X-uu ) 

at root of unity to = e 2m ' N intertwines not only the six- vertex L-operator, but also the 
following L-operators containing Weyl elements X, Z (see also |U]): 

/ 1 + AfoxZ; \X-\ ai -b 2 Z) \ 
L(X, a) = ; XZ = cuZX, (21) 

\X(a 2 -6 3 Z); Aa x a 2 + b 2 b- i b^ 1 Z J 

where A, a±, . . . , 63 G C . We collectively denote the parameters ai, . . . , 63 by a. In 
the representation © the intertwining relation is 



^ lil> ^ 2 (A^)L^(A, a )Lgf 2 (^a) = ^(^a)LS(A,a)^ lfaj2fc2 (A,i/).(22) 
iu'2,/3 ji j'2,/3 

where greek indices run over the values 0, 1, . . . , N — 1 and the latin indices take the 
values 0, 1 , #11,22 = w(A — z/), #12,21 = A(l — u), etc. 

Moreover, Bazhanov and Stroganov found that there is also an intertwining relation 
with respect to the Z N (greek) indices, i.e. in the Weyl quantum space, if the parameters 
a are chosen as 

viVi' . l vi' . u y<i 



a\ = x q ; a 2 = ux q i] b\ 



(23) 



HqfJ-qi ^q^q' VqVq' 

where the x q , y q , \i q etc. satisfy the CP conditions (jTSjl with fixed k. Writing (|2~T|) with 
the parameters (|2~5j) as L(A; q, q') , we get 



/ 



L(X; q,q') 



1 + A 



X I ClJXg/ 



Vq , 



AX- 1 

A US XnXgl + 



q^q' ) 



1 



(24) 



/ 



Apart from the spectral parameter A, this L(A; g, q') depends on three independent 
continous variables, e.g. x q , x q > and the modulus k . We shall not write the latter 
explicitly as an argument. The quantum space intertwining relation is 

E SwJk*M>Q><fi Li^i^P') L^ 2 (X;q,q') 

= E L Sf (A;?,9')^f (A;P,P0 S ft)fe7l)72 (p,p', g ,g')- (25) 

The matrix S turns out to be the product of four CP-Boltzmann weights (j!7j) : 
S aia2 ,p lf 3 2 {p,p',q,q') = W pq ,{ai- a 2 )W p > q {l3 2 - j3i)W pq {l3 2 - a^Wp'q'iPi- a 2 ). (26) 
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One can verify the relations (|25j). by explicit calculations, using (0) and (fT7)l several 
times, e.g. 

^ift b /3i./32;7i,72 - b Qi-l,/32;7i I 72 ~ A^/V ^ _ UJ o n -a2 Xp y q - ujfo-m+ly^ b ai,/32;7i,72- 

The Bazhanov-Stroganov periodic quantum chain of length Q is defined via its 
L-operator L^f(A; q,q') given in (|24jl and the corresponding monodromy matrix 

M (A, {ft, ft}^ 1 ) = L {^Q^%)L{\]qi,q' l )L(\]q 2 ,q' 2 ) ... L(X; q Q _ u q^), (27) 

where each L has its pair of rapidities qi,q[ and all these rapidities may be different 
while keeping the Baxter modulus k to be the same for all L. The transfer matrix is 

t = Tr C 2 M . (28) 

Baxter [TOJ ED] calls this model the T2(t g )-model after the notation r 2 introduced for 
the transfer matrix in equation (5.33) of 8J. Mostly in Baxter's work not the fully 
inhomogenous model is used, but the rapidities take two alternating values. In (jA.7|) we 
give the relation to Baxter's notation. Fusion of the BS-transfer matrices leads to the 
functional relations mentioned in the last subsection. 



2. 3D approach to the BS model 

2.1. L-operator 

We now show how the BS-L-operator can be obtained from the Linear Problem Q of 
the 3D approach, imposing a periodicity condition. We follow the general argument of 
|24j . which also gives a quantum group background of the construction. 

Consider the domain of the auxiliary plane containing the four variables 
roi, rDi, ro 2 , tt) 2 , see Fig. El The co-currents around each Weyl variable are taken 
to be related by the Linear Problem (|4"j). We impose the periodicity condition 

<V-i| = e^il; <0_i|=£(0il; <x-il = e<xil (29) 

in the vertical direction, with £ a quasi- momentum. Then the conditions (jlj) at fDi and 
rDi (those in the left hand dotted box of Fig. El denoted L\) are 

= < -0o| +£^ 1/2 (^i|ui + (0 o |wi +f«a(0i|uiwi , 

V ) 

= UJ 1/2 ( V->o| u l + ( 0l|Wi + Ki(0 o | u l w l- 

These linear relations can be rewritten in matrix form as follows 

(VIKmui - i)wr x = < 0| (3i) 

where (</>| and {ip\ denote the rows ((0o|> (4>i\) an d ((V'qIj ("01 1); respectively, and Ia(£) 
is the following 2x2 matrix with operator valued elements 

^ 1 — w^^UiiiiKiWiW^ 1 — Ui (a// 2 — KiWiW^f 1 ) ' 



MO 



l £ui («i — w'^WiWj x ) — u; 1//2 £uiUifi;i + wiw 1 1 i 



.(32) 



Bazhanov-Stroganov model from 3D approach 



9 







y 


e<xil 


<V>o| 







(Xol 

tt) 2 


-c 

M 

i 


^1 


1 

1 


(Xil 

^2 



Figure 3. Piece of the auxiliary plane which corresponds to the product of two L- 
operators. The Weyl pairs together with parameters K\, «2, Ki, £2 are associated with 
the corresponding vertices in this plane. 



We want to use a matrix representation for the Weyl elements. Observing that only the 
three elements wxwf 1 , u 1; ux appear, we may use a iV-dimensional representation, 
writing 

u x = u\ X ; ux = u\ X -1 (33) 



WxWx 



in the basis In order to bring this into a form comparable with (|24|) we put 

_ l/2 y^_ Wl _ _x y q y q > ^ 



K\ = UJ 



1/2 x q 

y q > 



so that 



( 



1 — ^U\U\- 



Ki = 00 



Ik 



X q ' 



Wl _ 1 VqVq' 
Wl XqXq'jJLqjJLq* 



Xq'/JgHqi 



-MxX ( UJ 



,1/2 _ y_3_ 



\ 



-1/2 Vq -^-1 



-E,uiUi h 



Xq'Hq^q' 
VqVq' 



\ 



(35) 



X g ' OJXgXq'flqUq' ) 



XqUqUqi 

Conjugating with the Pauli matrix o~i and introducing a new spectral parameter A by 

1 



A 



we obtain 



oa 0-2 



/ 



A UJ XqX q i 



1 + A^Z: 



X 



Aa; 1,/2 tiiX„X 



V 



■Ux£„A. a; x 



flqflq' 



Vq 



fXqfXq' 



UJ 1 l 2 U\X q 



A a; + 



/iq/V 
1 „ 



A gauge transformation with 

, \/u) X l 2 UxXq\ 







1/ \fuj 1 l 2 uix q \ 



(36) 
\ 
/ 

(37) 
(38) 
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leads to 

P x a 2 L x (0 °2 Pf 1 = — ^ L(A; g, g') (39) 

\UJX q y q 

with L(A; g, g') defined in (J23j) . We shall see in (j3Hj) that the identification (J53j) will also 
appear when we express the BS intertwining matrix S within the 3D framework. Using 
the parametrization introduced later in section 1231 we will show in (|7U|) that considering 
the monodromy matrix, the factor U\U\X q y q relating the spectral parameters £ and 
A does not depend on the site considered. The same holds for the combination u\x q 
appearing in (j38j) . So the BS-monodromy (J27j) can be written as 

Q-l 

M (xAQu&to 1 ) = Px^MOHO ... Lg^^Pf 1 \{ (Xojx qi y qi ). (40) 

i=0 

We finally remark that demanding periodicity (ffilj) not after two vertical steps as 
done here, but after iV steps, one obtains N x N L-matrices, see j24j . 

2.2. 3D interpretation of the relation SLL = LL S . 

Consider the product of the successive action two L-operators (the simplest monodromy) 

<?/>|(^u 2 u 2 -lH^mCu -ljw^wr 1 =( x \-L 2 (OL 1 (0. (41) 

We are interested in the relation of the action of L 2 {^) -£a(£) to the action of 
Li(£) L 2 (£). In the 3D approach, 1Zi 2 % maps a triangle into a reflected one, recall the 
left hand picture of Fig. [T] Let us consider the four Weyl operators in the auxiliary 
plane as in Fig. El (the scalar Kj is included in tVj, see 0) and a further variable W3, 
see the bottom auxiliary plane in Fig. |3J 

We consider the special case of 7Z\ 2 3 mapping the triple Weyl algebra ( tt>i, K> 2 , tn 3 ) 
into (tt>i, tt> 2 , tn 3 ) with the functional part of the transformation TZ 123 taken to be 
trivial (Bazhanov-Baxter case): v! i = Ui, w[ — Wi, i = 1,2,3. We act with a similar 
mapping TZ-q 3 on the initial triple Weyl algebra (fOi, tt> 2 , tv' 3 ) and obtain the triple 
algebra (ro^, ro 2 , tUg') . As in Fig. |3]we demand periodicity after the second step in the 
third direction: tog = tt> 3 . Fig. |U illustrates these mappings. 

In the upper auxiliary plane the action of the operators L\ and L 2 appears in the 
reversed order if we keep track of the direction of the lines. 

So, taking into account the property (|SJ) we expect an intertwining relation of the 
form (|25p . with S being a bilinear expression in R123 and Rygg. 

2. 3. BS matrix S from ZBB model 

In order to derive the precise relation, let us calculate the trace of the product of two 
matrices R and R, both of the form (jllj) . but the first with the Fermat parameters 
Pi = (xi,yi), the second with pi = (xi,yi) , i — 1,2,3, 4, each satisfying the restriction 
©: 

ciii2,fcife 2 _ \ ^ -piii2« -Bfcifem (ao\ 
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Figure 4. Graphical image of the intertwining relation for the BS model showing the 
origin of the Bazhanov-Stroganov intertwining matrix S. The two elements roi and 
fbi form the operator la (A), the elements tt>2 and tt>2 the operator L2(A), see the 
earlier Fig. [3] Periodicity in the third direction gives that the two tt>3 at the top right 
and bottom left are the same. The two points where IZ123 and ac t are vertices 

of the physical 3-dim lattice. 

Actually, we shall not need the full iV 4 x iV 4 matrix (J4*2|) but only those matrix elements 
in (g2J) with 

ii + £1 = %2 + h = Ji + h = j 2 + k 2 = 0, (43) 

forming a N 2 x N 2 matrix. Inserting the explicit expressions ()11|) and renaming the 
discrete variables we obtain 

S Pi02 = w (oi-/3i)(A-a a ) w Pl (a 2 -a^wpjjxi - a 2 ) Wpjjh - fa)wpjjh - (3 2 ) 

aia2 " w Pz{(^2 - ai)wp 3 (o;i - (5 2 ) w P4 (a 2 - (3i)wp 4 (f3i - a 2 ) ' 

In order to rewrite (J44j) in the form (|26|) we use the following property of the cyclic 
functions w p {n)\: 

w p( n ) = / ~ \ /r, / \ , (45) 

where O is an automorphism of the Fermat curve such that 

p = (x,y) 1— > Op = [u^x" 1 , uj^^x^y) (46) 
\ Relation (|45|) is an analog of well known relation r(x)r(l — x) = — , see e.g. |4l 113). 
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and 



$( n ) = (-) n oo n2/2 . (47) 



So we get 



^ Wp^ai - a 2 ) w P2 (/^2 - A) w p 3 (f3 2 - «i) WQpJJh - a 2 ) 

w 0pi {a£i - a 2 ) wopzifo - Pi) w P3 {f3 2 -ai) Wp 4 (p 1 -a 2 ) 
Now we can identify 

u^Oi-c^) TJr ^ w P2 ((3 2 -(3x) 



(48) 



wop! (ai - «2j w p 2 (P2-pi) 
vv pq (p 2 -ai) = t- -; H/p/ g /(pi - a 2 J 



(49) 



w P3 (/? 2 -ai) ' " " WpJJk-012) 

because of the trivial identity 

rcK-Mifr-az) $(«l-/3 2 )$(« 2 -/3i) = 

$(ai-a 2 )$CS2-A) ' 
The notation for the rapidities p,q,p',q' which parameterize the point on the Baxter 
curve (J18|) should not be mixed up with that of the eight Fermat curve points 
Pi, pi (i = 1, 2, 3, 4) . Using flUJ) and (JT7J), (HH we can write 



www 


JW 1 ■ 


-u n {x p /y ql ) _ 


uj x l 2 y\X\ 1 — uj ti 


(wxi) 1 


W p , q ,(n-l) 


Vq'Vp 1 ■ 


- uj n (x ql /y p ) 


yi i- 




W p',q( n ) 




-uj n (x p/ /y q ) _ 


(jj 1/2 y 2 x 2 l-uj n 




W P ', q (n - 1) 


/W 1 


-uj n (x q /y p/ ) 


m 1 - 




and 










W M (n) 




1 - u n ~ l (x q /xp) 


_ ^" 1/2 y 3 1 


— u; n ~ 1 (t<;x3) 


W M (n ~ 1) 


x p 1 Uq 


1 - uj n (y p /y q ) 


2/3^3 1 


— ^(uxz)^ 1 


W p ' jq '(n) 


UjXpflXqi 


1 — uj n ~ x {x q i 1 Xp 


_ ^" 1/2 Z/4 1 


— u; n_1 (u;x4) 


W p ' )q '(n - 1) 


X ~phq> 


l-uj n (y p i/y q/ ) 


2/4^4 1 


— 0J n {tjJX4)~ l 



(51) 



(52) 



From (|51|) and (|52j) we can read off the relation of the Fermat parameters to the CP- 
variables on the Baxter curve: 

Vq' X q %q Vq' 

Xi = ; x 2 = — ; x 3 = —2- ; x 4 - 



uXp y p i ujXp ujy p i 

(53) 

x q' ~ Uq ~ Uq x q' 

X\ — — ; x 2 — 5 X3 — ; X4 



y p UXp/ uy p ujx p i 

Note that both restrictions (fT2|) are valid automatically for the Fermat points pi and 
pi. Comparing the coefficients in front of the ratios in (joTj) and ()52j) we can identify the 
ratios fji/yi with ratios of the Baxter curve parameters 
yi = 1/2 [ipXp _ ^2 _ 1/2 /W . 

2/1 /y%> ' 2/2 /yay ' 

§3 = w l/ 2 /Wp . ?/4 _ „-l/2 2/p' 
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From (J54j) we obtain an important relation which connects the Fermat points pi and pf. 

ym mm 

= ■ (55) 

2/31/4 2/32/4 

By taking the iV-th powers of the formulas f)54j) one can see that the cyclic property 
of the functions w p (n) (jl(Jj) implies the cyclic property (J 19)) of the Boltzmann weights 

(ED 

We conclude this discussion of the emergence of the BS-S-matrix within the ZBB- 
model with showing that the parametrization of the L-operator postulated in (}34j) 
agrees with the identifications (|53j) . (|54|) made here. First we pass from the parameters 
Ui, Ui, Wi, Wi, Ki, hi used in to the Fermat parameters by (JT3J). Since the functional 
transformation is taken to be trivial, we can omit all primes on the Ui and u>j. 

-1/2 u 2 -1/2^2^2 _l^2 

x\ = uj 1 ; x 2 = uj ' ; x% = uj — ; 

2/3 Kl^l 2/4 -1/2^3^3 ( r R \ 

— = ; — = uj 1 . (56) 

2/i «3 2/i w 2 

The counterpart with tildes is: 

_l/2 U 2 _ _i/ 2 /«2M2 ~ _l«2 

Xi — uj 1 - - ; x 2 = uj 1 — — ; x% = uj — ; 

K\U\ U\ Ui 

2/3 KlWl 2/4 -1/2^3^3 , K -v 

— = ; — = uj ' — . (57) 

2/i u 3 y 1 w 2 

Observe that there are no tildes on k 3 , 113 and 103. (|56[) and ()57|) immediately give 

1/2 X 3 ~ 1/2^3 Wi Kl 2/1 2/3 

xi xi wi «i yi y 3 

M 2 u 2 „ 

— = 07x3 ; — = 1^X3. (58) 

The first three equations of (|3*4"|) follow by simply inserting from into the first three 
equations of (|5H|). The last two equations of with (|S3*j) lead to 

x p u 2 = x q ui and y p u 2 = y q u\ or Uxu x x q y q = u 2 u 2 x p y p . (59) 

fj59j) shows that the rescaling of the spectral parameter = —uju\U\x q y q A is the 
same for L 2 (X;p,p') as it is for Li(X;q,q'). 

Summarizing, the relation of the parameters K\, K\, k 2 , k 2 , —, — to the CP- 

Wi w 2 

parameters is: 

1/2 x q ~ -1/2 2/g i/2 x p ~ -1/2 2/p 

k x = uj 1 — ; K\ = uj ' — -; K2=uj' — ; /t 2 = c<j 7 — -; 

2/g' ^g' 2/p' V 

^ = W - 1 _JW 2 L_ !^ = fa ;-i ^ . (60) 

U)l XgXg'fiqfiq' W 2 X p X p >fl p (i p > 

2. 4- Parametrization in terms of cross-ratios 



The intertwining matrix S defined by (j^fij) depends on 5 independent continuous 
parameters. One may use several equivalent parameterizations: 

(i) CP-parametrization: q, q', p, p', k . 
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(ii) Fermat parametrization: Xx, x 2 , x 3 , Xi, x 2 , x 3 with the constraint (|55|). 

(iii) Weyl-parametrization: k%, Kx, k 2 , k 2 , Wx/wx, w 2 /w 2 , again with one constraint. 

One obtains (ii) from (i) by (J53J) and (|54|) . while (iii) follows from (i) by ([60)1 . In (iii) 
one may choose u 2 /ux , u 2 /ui instead of wx/wx, w 2 /w 2 as is seen from (}56|) . (j57j) . 

In the following let us concentrate on the functional transformations, which 
deal with only the N-th powers of the variables. We start considering the Weyl 
parametrization (iii). The functional transformations have been solved in j2Z|, first 
rewriting them in Hirota form and then using standard methods of algebraic geometry. 
The results can be read off from (65), (66) of [23] and Table 1 of [2Z1- m the moment 
we are dealing with only the trivial version of the functional transformations. Their 
solutions can be obtained either by specializing Table 1 of |27j . or directly from (j!4|) 
solving the system 

{K 2 U 2 - KJJ X ) W 2 = {K X U X + U 2 ) K 3 W 3 ; 

(Wx - K 3 U 3 ) W 3 = (Wx + U 3 ) W 2 ; (Ux - U 2 ) U 3 = {K x Ux + U 2 ) W x (61) 

in terms of three pairs of complex points which we shall call X', X, Y', Y; Z' , Z . 
The solution is 

I; u 2 =-e^—^- u 3 



Y -Z, 







Wx = e\- W 2 = W, 



Y-Z 







where 

e = (-1) N ; Kt = «f, i = l,2,3, (63) 

and for cross-ratios we use the notation 




(A-C) (B-D) 
(A-D) (B-C) • 



(64) 



There are analogous equations for the variables with tilde. We solve these introducing 
another pair of complex points Z[, Z±. The expressions are the same as (J62j) . only 
everywhere Z' , Z replaced by Z[, Z x , e.g. Ux = —e(Y — Z[)/(Y — Zx), etc. 
Observe that the ratios appearing in (iii) can be written as cross-ratios: 




(65) 



Wx I Y' Y I W2 

Wx \z Q z x y w 2 

so that all six variables in (iii) are expressible as cross-ratios of eight complex points. 
Frequently, we shall use calligrafic letters to denote pairs of points: 

X = (X',X); y={Y',Y); Z = (Z' ,Z ), Z x = (Z[,Z 1 ). (66) 
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Due to projective invariance, from these eight points five independent cross ratios can 
be formed. Observe the index 1 is associated to the pairs y, Z, the index 2 to the pairs 
Z, X and index 3 to X, y. 

The Fermat parametrization (ii) is obtained from (|62|) by (|56|h ()57|): 

n _ ) ^ 1 i . \ i 1 ^ i . v ; y x I 7v I x' y 

Xi — ^ „, „ ^ i Xn — < „ „, > , — < „ „, > , X 





and analogous formulas for Xj and yj where the pair Z is replaced by Z\. Observe that 
the Fermat curve conditions (fTU|) and (|T2")l are trivially satisfied. 

We shall see now that also Baxter's modules k 2 , k' 2 and the iV-th powers of the 
CP parameters x p , y p , etc. have similar good expressions in terms of points X, X', . . . 
Since (153)1 and (J5l|) involve only ratios of the CP variables, to solve e.g. for x^ and 
Up we have to use the Baxter curve relation (JTSJ). 
One may proceed as follows: From (JBOj) and (|52)l get 

~ ~ \ N 

K 2 W 2 \ , n TV 



K 2 W 2 




Then use (fT8|) to write 



Substitute here from (|50|l ay = a; 1 ' 2 y p / k 2 and solve for y. 



V^i //f.r. snrl snlvp for 11^ 



i . N _ 1 ~ (/V/v) 



1 + (/ip/v) X, 



2 



and insert, using also (|52)l . To obtain x^ the same procedure works. So we find all 
CP-variables: 



k Xp — \ ry rjf f ) ^ ~~ ^ ry r^l \ j ^ Z 1 ^ 




z x z; 



/ N _ J yy fl 1 t, N _ ) X' Z{ i 7 / jy 

fta V - 1 7, 7/ r ~ 1 zJq z^ ' " /,,/ 

Zi Z( 



1 N _ ) 1 ^0 I I. JV _ J ^ ^1 I i/ iV 

Kx i' ~ \ v. v< r > ~ i z z^ • ' ^ 




(6* 
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Each line of (jUH|) yields the same expression for k 2 = kx^ + ky^ — k 2 x^ y N 



k' 2 = < Z z [ \ \ or k' 2 = { 7 ? > . (69) 

We also note that Baxter's rapidities p, p', q, q' correspond directly to the points X, 
X', Y, Y' , respectively, while Baxter's module k does not depend on these points. 

Returning to the parametrization of L-operator (JH2I), in (|3*3*|) we found that it 
depends on three parameters K\, k\ and W\/w\. These are cross-ratios of six points 
Y, Y', Zq, Z' , Z\ and Z[. Due to projective invariance six points give rise to only three 
independent cross-ratios. These correspond to the rapidities q, q' and the module k in 
the Bazhanov-Stroganov parametrization (j24j) . 

Finally note that in this parametrization the combinations encountered in ()36j) and 
depend only on Z and Z 1 : 




^<< = f-rf; u '< = - e k§r I h- (70) 



3. Classical BS-model and intertwining of its L-operators 

3.1. Functional mapping on N -powers of Weyl operators 

As has been mentioned before, for u N = 1 the A^th powers of the Weyl operators are 
central and can be considered classical variables. As in (J7J) we use {i = 1, 2, 3) : 

E/i = uf; Wi = wf; U { = uf ; W t = wf ; K t = «f ; = «f , (71) 

and A = ^ 7V . 

According to the approach developed in |2E] we introduce a classical analog of 
the linear problem (J3*U|) 

= # - x A#i + $0^1 + SjA^itWi 

(72) 

= _ + + ^o^C/xWi 

which can be rewritten in the matrix form 

(Atfitfi - 1) = $ ■ £i(A) (73) 
and defines the classical L-operator 

( Wt + AU^KiWi Ux (Wx + KxW^j ^ 

v AUx (Wx + KxWx^j Wx + AKxUxUxWx > 

acting in the space of the linear variables \I> = (\I/o> ^i) and $ = (^o^i)- We take 
f!74|) to define a discrete classical analog of the Bazhanov-Stroganov model. We can also 
define (174)) by an averaging prescription from (|32j) : Define ^7] 



A(A) 



(74) 



i G Zjv 
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Then 

A(A) = W x ■ (L x (0). (76) 

Analogously, we introduce an operator £2 (A) such that the classical variables and 
parameters U\, U\, W\, W\, A 1; K\ are replaced by U 2 , U 2 , W 2 , W 2 , K 2 , K 2 . Let 
£*(A) and £ 2 (A) again be L-operators of the form (JUJ), but with the variables 

U*, W*, Uf, W?, Ki, Ki, i = l,2. (77) 

Our aim is now to find the transformation (Ui, Ui, W\, . . . , W 2 ) 1— > (U*, U*, W*, . . . , W 2 ) 
which solves the intertwining relation 

£ 2 (A)A(A) = £*(A)£*(A). (78) 

However, trying to find the nonlinear 8- variable mapping by direct calculation without 
further guidance looks quite hopeless. Fortunately, the 3D approach will provide a 
solution to this problem 

3.2. Solving the classical BS-intertwining relation via the 3D functional transformation 

The mapping (|78jl we are looking for can be found using the functional mapping of the 

vertex ZBB-model given in (JTlJl 5|. We introduce two additional variables U 3 , W3 and 

( f) 

the additional parameter A3 and consider the rational mapping lZy 23 

K[& : U x ,W t ,U 2 ,W 2 ,U 3 ,W 3 h- U[,W[,U' 2 ,W^,U' 3 ,W^ (79) 

given explicitly in t)14j) . We define the composition of two of these rational 
transformations f|79jl 

n { / 2 \ ■. u x ,w x ,u 2 ,w 2 ,u z ,w z ut,wt,u*,w*,u' z ,wi, (80) 

74§ 3 : f/i,^,^,^,^,^ h- ^,T^,E/£,W?,E/* (81) 
together with a periodic condition 

U£ = U 3 , W* = W 3 (82) 
and denote this composition by 

S$ : u u w 1} u 2i w 2i u u w 1} u^w 2 

^ Ut,W*,U*,W*,UlW*,U*,W* . (83) 

In ()8L)|) the constants Ki,K 2 , A 3 have to be used, while in (J81)) the constants are A 1; K 2 
and A3. With these definitions we have 

Proposition 1 The rational transformation S\ 2 \83\) solves the mapping defined by the 
intertwining relations \78j) . 
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The Proof is provided by straightforward calculation. We first determine U£, W£ in 
terms of the variables U 3 , W 3 and other variables from the successive application of first 
(|8T|l and then (fHUJ) . Imposing the periodicity condition (fH2~|) gives two equations which 
can be solved easily for the auxiliary variables U 3 and W 3 , leading to 

U 1 (K 1 U 1 W 1 + K 1 U 2 W 1 ) + U 2 (U 2 W 1 + U X W X ) 



U X U X - U 2 U 2 
W 2 W 2 {K 2 K 2 U 2 U 2 - KxKxUxUx) 



(84) 



Then 



K Z (XJ 2 (K X U X W 2 + K 2 U 2 W 2 ) + K X U X (K X U X W 2 + K 2 U 2 W 2 )) 
is used to eliminate U 3 and W 3 and we get 



Ut 



u* 



TT* 



k x \J x {K x U x W 2 + K 2 U 2 W 2 ) + ggCgigi W 2 + K 2 U 2 W 2 ) 

' K 1 u 1 (k 1 u 1 w 2 + k 2 u 2 w 2 ) + u 2 {k 1 u 1 w 2 + k 2 u 2 w 2 ) 

ih Utik&Wi + k x u 2 w x ) + u 2 {u 2 w x + u 1 w 1 ) _ 
us ~ 



w x w 2 v 

where 

V = {K x KxUxUx - K 2 K 2 U 2 U 2 ) {UxU x - U 2 U 2 



U 2 W X (K X U X + U 2 ) + UxWxiKxUx + U 2 ) 

K 3 Vx W[_ _ W 2 K 3 

Wx 



V 



WxW 2 V 



55) 



(86) 



Vx = kxUxlll [U 2 {Wx + KxWx){W 2 + K 2 W 2 ) + K x U x WxW 2 

+ K 2 U 2 U\ \Ux(Wx + KxWx){W 2 + K 2 W 2 ) + K 2 U 2 W X W 2 

+ UxU x U 2 U 2 W 2 W X W 2 {\ + KxKx) + KxWxW 2 {\ + K 2 K 2 ) + 2 K x K 2 WxW 2 

and V2 is obtained from Vx interchanging the variables with tildes and without tildes. 

Finally, we simply insert into (|78*|) . These are 2x2 matrices with four entries and we 
compare the coefficients of A , A 1 , A 2 each, giving 12 equations which turn out to be 
correct. □ 

Because of the periodic boundary conditions (|K2*|) the mapping (|K3|) has the invariants 

UxUx = U^Ut ; U 2 U 2 = U 2 U 2 ; W X W 2 = W{W 2 * ; W X W 2 = W{W 2 . (87) 

Note that the first two invariants in ()87)1 reflect the fact that the combinations UU are 
scale factors of the spectral parameter, see (JoTJj) . The last two invariants show that a 
difference in the normalization of the classical and quantum L-operators (fT4"j) and (}32|) 
is not important. 



4. Main result: Isospectral transform of the BS transfer matrix. 

We now consider the BS-quantum chain of length Q defined by the monodromy 

M (0 = L o(^ u , . . . , k ) , ux,...,Kx) ... Lq-x(£, uq-x, • • • , kq-x) 
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and the transfer matrix 

T(0 = tr c2 M(f) 
where the Lax operators are defined by (JH2J), 



/l 1/2 1 ~ n ry V f 1/2 ffl H ry 

/ 1 - w 7 £u n u n K n — Z n ; -m„X„ u; 7 - « n — Z n 



9) 



^ i / 1/2 ^^'ft> \ 1/2 ~ ~ 

i £ U n X„ \K n - LU 1 —Z n ); -U) 1 £u n U n K n + — Z n/ 
\ V Wn J W n / 

The arguments tt n , «„, w n /w n , n n , k n of all P n (n = 0, . . . , Q — 1) may be different 
and shall be parameterized following (ffi2*|) in terms of Q + 2 pairs of points y n and 
Zo, Z\ : 

In the equivalent formulation of the transfer matrix in (|2*H|) in terms of CP variables 
x qo , x q ' o , . . . , k 2 this means that we take the Baxter modulus k to be the same for all 
L n but the two rapidities in each L n may be different. The normalization adopted in 
dHHJ) differs from (|27|) . (|2gjl . see 

The main result of this paper, which will be proven in the remaining part of this 
section, is the following: 

Proposition 2 For a given Baxter modulus k and a fixed set of 2Q rapidities, there 
exists a Q — 1 parametric family of transfer matrices with the same spectrum as 
the initial one. This family of transfer matrices is defined by the same formulas 
$EB), (GSJ where the matrices satisfy as usual X ni Z„ 2 = o/ n i n 2 Z„ 2 X ni , but 
without being normalized to = Z„ = 1 . The form of the inhomogeneous centers 
is given by Hlty) . 

In order to state the results for X^ , Z^ , we have to consider the following system of 
algebraic equations for the unknowns P, P' , with the pairs 3^n; ^i being given: 

Q-l 

II I-' !-' I 

1 ; (91) 



1 . (92) 




i=0,l 



The system (f9T|l . (J92|) has exactly g = Q — 1 nonequivalent solutions {(P-, P-)}, 
j = 0, g — 1, the pair (Pj, Pj) taken to be equivalent to the pair (P-, Pj). 

Now, given a fixed set of g pairs of complex numbers (Pq,P ), . . . , (P' 1} P g -x) , 
we define the function H of a g- dimensional vector (/ , fi, f g -i) = {fj} , denoted 
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H ({fj}) (this arises as the rational limit of the B-function on a genus g generic algebraic 
curve, see the Appendix of [15]) § by 

Apt I pi — f.P'i 

'«*»- kMr ' (93) 

This is seen to be normalized to H({0}) = 1. Let further 

UY) = ffi a s (y) = *,{¥',¥) = j J }> (94) 

recall © 3> = (V, ^) , and define 

iM) = H v vW = II ^O*); W) = i s {Q) = i- (95) 

i=o I * » J i=0 
Then, the spectrum of the transfer matrix defined by (|88|h ()89|) . (j9~Uj) with 



H({^(Y n )I ] {n)a 3 (Z Q )}) ' 
# ( mZ ) I 3 {n + 1)} ) ( ft(Z0 J>(n)} 



(96) 



( {^(ZO J i (n+ 1)} ) ( {f,-(Zo) /.(n)} ) 



does not depend on the set {fj}. If {fj} = {0}, then (|9l)Jl reduces to = = 1 , 
and so the spectrum is the same as for the initial transfer matrix. 

In order to prove these statements we shall make extensive use of the 3D-formalism 
which will allow an easy description of the necessary intertwining operations. We shall 
derive a generalization of the parametrization (J9*U|) which will involve the functions H 
which depend on the parameters {fj} . A key point will be that the product u n u n 
and the K n , k n will not involve the {fj} , see eqs. (jl()8)) . and that in ()88jl Z n is always 
multiplied by (the {/jj-dependent) factor w n /w n and X n is multiplied by the {/■,}- 
dependent u n . Details should become clear as we proceed. 



4-1. Uniformization of the classical maps 

The map ()83j) describes the explicit relation between the final "star" variables and the 
initial "non-star" variables. According to [21 121] this map, as well as map (JHJ), can be 
parameterized in terms of algebraic geometry data. In this paper we will consider the 
uniformization of these maps using a specific set of the rational functions and identities 
between them. This construction was exploited previously in |15j . 

Consider a general three-dimensional lattice with classical variables placed on its 
edges. Let n = niei + n 2 e 2 + n 3 e 3 be marks for the vertices in this 3D lattice with 

§ These rational functions appear in soliton theory (see, for example, |14| . g: number of solitons), and 
in 3D integrable models |3()j : 

H( ) = 1; H(f ) = 1 - /o; H(f , h) = 1 - - + ^ foff, 

H(f f f \ - 1 (Pj-P^(Pj-Pi) f (P{-P 2 )(PZ-P2){Pl-PZ) f f , (P;,-P[)(Pi-Pi)(Pi-Pa) f f f 

nU0,JUj2) - (p _ Pl )(p _p 2 ) Jo^- ■ •+(P 2 -P 1 )(P 2 -P )(P 1 -P ) 70/1-I---- (P 2 -P 1 )(P 1 -P a )(P 1 -P ) /0/1/2- 
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U 2 ,n, W 2 



Ui, n , W 1>a 



"l,n+ei; Wl,n+ei 




2,n+e 2 ? »' ,/ 2,n+e2 



u 3 , n , w 3 „ 



Figure 5. Association of the dynamical variables to the edges of the cubic lattice. 
The vector ei is pointing into the paper plane. 



discrete coordinates rii,n2,n 3 . The assignment of the classical variables to the links 
around a given vertex n is shown in Fig. |5J In this notation the map (|14|) relating the 
neighboring classical variables is 
Ui >n+ei W 3 



3,n 



K 2 :ni,n 3 U 2 , a W 2 ,n 



:ri2,n 3 U ltn W 2tn + K 3:ni>n ,U 2 , a W 3>ri + Kt K 3 

Win W 2 , n+e2 W hn W 3 „ ( 9? ) 



Wl.n+ei 
^2,n+e 2 



'2,n+e 2 

w 2 , n 
u 3 , n 



W 1>n W 2>n + U 3 , n W 2 , n + K 3 

:ni,n 2 

u 3 , n w 3!n 

u ljU u 3jn 



u 2 , a u 3 U 2 ,JJ 3 , a + U 2 , n W 1<n + K x 

One may think about these relations as discrete equations which describe the 
interrelation of the classical variables along the 3D lattice. 

The next step is to observe that after the change of variables 



u. 



2,n 



U; 



3,n 



— e 



Y 

1 n 2 


- Z 'n 3 


T 2 ,n 


Y 

1 n 2 


— z nz 


T 2 ,n+e :i 


x m 




Tl, n 


X m 


— Z n . A 




X m 


-Y' 
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T i,n+e 2 


X ni 


_ Y 

1 n 2 


Tl, n 



2.n 



W: 



3,n 



Zn 3 


^n 2 T 3,n+e 2 


Zn 3 




Zn 3 


- X' ni r 3) n 


Zn 3 


X ni T3,n+ei 


Y 

1 n 2 


-X' ni T 2 ,n 


Y 

1 n 2 


— X ni T~2,n+ei 



together with the cross-ratio parametrization of the Ki- n . nk 

ly l:n 2 ,n 3 \ ryi ry (i * l 2:ni,n 3 ) ^ ( j - fv 3:ni,n 2 " 



y/ y 

1 n 2 1 n 2 
Z'n 3 Z n . A 



7' 7 

n 3 n n 3 

x' ni x ni 



x' ni X 



"i 



Y> y 

1 n 2 1 ri 2 



(98) 



, (99) 
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each of the relations in (|97|) can be written in the form of a three-linear Hirota-type 
equation for the triple of unknown functions r a]I1 , a = 1,2,3: 

(x a — Xfs)(x' /3 — x! y )(x y — x a ) 

+ (X a -X' )(X -X J )(X! y -X 

= (X a — Xp)(Xj 3 — X 1 ){X f 1 — Xajr^n^r^n^r^n 

+ (X a -X' fj )(X -X! y )(X J -X 

where {a,P,^} is any even permutation of the set {1,2,3}. The notations in ()100|) are 
related to those in (J98|) as follows: 

X 1 = X ni , X 2 = y' n2 , X 3 = Z n3J where y' = (Y,Y') ■ (101) 

Now our goal is to describe the general solution to (jlOOj) in the ring of the "rational 
G-functions" H defined in (|93|). The main tool to be used will be an identity, which is 
the rational limit of the Fay identity ("rational Fay-identity"): 

Let A, B, C, D be any pair-wise different complex parameters. Using the definitions 
given in (JHBJ),(jniJ), the following identity is valid: 




H({f 3 }) H({f 3 a 3 (A } B)a 3 (C,D)}) = { n n } H ({/^(A, B)}) H ({f^C, D)}) 




+ < u n \H{{f 3 a 3 {A,D)}) H ({f 3 a 3 (C, B)}) (102) 



For a proof of this identity see the Appendix to the paper ^H] ■ In order to get the form 
of ([IDOp . we combine two such identities and obtain the more complicated "rational 
double-Fay identity" 

(X-Y)(Y'-Z')(Z-X)H({f j (X)a j (y)a j (Z)}) H ({f 3 (Y)}) H ({f 3 (Z)}) 

+ (X-Y')(Y-Z)(Z'-X)H({f j (X)}) H ({fj(Y)o-j(Z)}) H ({f 3 (Z)a 3 (y)}) 

= {X-Y){Y'-Z){Z'-X)H{{i j {X)a j {y)}) H {{i 3 {Y)a 3 {Z)}) F({f,(Z)}) 

+ (X-Y')(Y-Z')(Z-X)H({f j (Z)a j (y)}) H ({f 3 (Y)}) H ({f 3 (X)a 3 {Z)}). (103) 

This identity involves five complex points X, y, Z and the sets {f 3 } and {Pj, P'}. 
Dividing by (X - Z'){Y - Z)(X - Y')(Z' - Y)/(Y - Z') the factors multiplying the 
functions H can be written as cross-ratios. The structure of both 1)102)1 and (J103j) 
is precisely the same as that of the corresponding identities for 0-functions, see e.g. 
equations (21)-(24) of 

Comparing ()100|) and ()103|) . we conclude that for arbitrary {/,} the discrete 
equations (jl00|) on the 3D cubic lattice are solved by (we have to use y' instead of 
y in ()103|) which explains the inverse in the middle term of ()105|0 : 

r ljn = H ({f, (X ni ) o.iy, , ,)!,..„} ) ; r 3 , n = H ({f, (Z n3 ) a 3 (y n2 )I j:a }) , 

r 2 , n = H ({f- {Y'J aMly.n}) = H ({fj (Y n2 ) I j:n }) , (104) 
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Figure 6. Slice of the 3-dimensional lattice associated with the BS chain. The 
orientation of ei, e2, e3 is shown in the left lower corner. The lattice is taken periodic 
after two steps in the e3-direction. The classical variables are assigned to the links 
as indicated in Fig. In the Figure we show a few of these variables where we just 
indicate the Ui tI1 variables, the Wi^ n at the same links are not shown. Observe that 



where we introduced 

Tlx— 1 712 — 1 "3 — 1 

I jm = H a 3 {X mi ) J] aj\y m2 ) J] a,(Z m3 ); I j 6 = 1. (105) 

mi=0 7712=0 7713=0 

This is the rational analog to solving trilinear Hirota equations by use of the double Fay 
identity for O-functions [HI Ell ■ 

Inserting (jl04j) and (jl05j) into we get the parametrization for all variables on 
the lattice: 



U hn = - e 



U 2 , n = ~ e 



Y n2 -Z' n3 H (ft (Y n2 ) I jm }) 



Yn 2 — Zn a H ({fj (Y n2 ) Ij. n O-j(Z n3 )}) 

X ni - Z' n3 H ({f 3 - (X ni ) I r . n cr,(y n2 )}) 
X ni - Z n:i H ({fj (X ni ) Ij: n a- j {y n2 )o- j {Zn 3 )}) 



rj X ni -Y> 2 H (ft (X ni ) I j:n }) 

3 ' n " X ni - Y n2 H ({f, (X ni ) Ir. a aj(y n2 )}) " 1 > 

The Wi: n are obtained from the Ui :n interchanging the arguments and removing the 
overall minus signs. 

4-2. The classical BS chain 

In the last subsection we considered the whole 3D lattice. In the following, we shall 
be interested in applying this 3D formalism to the BS chain defined by the product of 
operators C n . The chain will be taken in the e 2 -direction, see Fig. |3 As in Fig. |SJ 
for each C n we have to consider a pair of vertices, here on top of each other in the 
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e 3 -direction, with periodic boundary conditions after two steps. Just one layer with 
open b.c. is needed in the ex-direction. 

To each line of the 3D lattice we associate two pairs of points as shown in Fig. |HJ a 
line in direction e x is labeled by the two pairs 3^, Zj. Classical variables are associated to 
the links as in Fig. EJ In the notation of (J73|) we have Ui iXl = Ui, n+e3 , Wi yXl = Wi jn + e3 . 
The indices take the values n x — 0, n-i — 0, 1, . . . , Q — 1, n 3 = 0, 1 starting from the 
left bottom. 

Looking to the allocation of the variables to the links of the lattice, the Lax 
operators associated to the lines y n2 will be taken to be of the form (|74^1 and depending 
on the index " 1" -variables U\^ n2e2} Wi. n2e2 , i^i,n 2 and their tilde counterparts. Only n 2 
will occur. We drop the index 2 and write C\^ n . Explicitly: 

^ Wi^ ne2 +AUi^ ne2 Ui tne2 Ki. n Wi jne2 U\.nei (w Xjne2 + Ki. n Wi, ne2 ^j ^ 

^ Af/i ;ne2 (Wi, ne2 + K\:riW\,ne2 ) W\^ ne2 + kKx, n U\, n eJj\,ne 2 W^ ne2 j 

Since these classical variables are solutions of (jl00|) . according to (jl06|) we can write 
them in terms of the rational O-functions H defined in (jHSJ): 

j j _ Yn-Z', Hi^^Wm fj (Y n -Z>)(Y n -Z[) 1 

1,ne2 ~~ Yn - Z H{{ij{Y n ) Ij(n)aj(Z )}) ' 1>ne2 " (Y n - Z )(Y n - Z x ) U x ^ 2 ' 

^i,ne 2 _ jY; Y n \H({^(Z )I 3 (n + l)}) H({fj(Z )Ij(n)}) . 



A,n(A) 



(107) 



W hne2 \Z Z x \ HaWZMin)}) HdUZjIjin + l)}) 



K - -< ] zi Y iy \z\ l>- « 

Here {fj} is an arbitrary set of Q complex parameters. The Ij(n) are given in (J95j) 
and are related to the Ij. n of ()105J) by Ij(n) = lj-ne 2 ■ 

We define the classical monodromy matrix, the classical counterpart of (jHHj) . as: 

M X (A) = (]i W^J £ 1)0 (A) A,i(A) • ■ • A,n(A) • • • A,q-i(A). (109) 

\n'=0 / 

The simplest choice, still compatible with the functional mapping relating neighboring 
variables, is to take all {fj} = {0}, resulting in all H being unity. Then M. x still depends 
on the Q pairs y n , the two pairs Zq, Z\ and the spectral parameter A. Since all y n may 
be chosen differently, in general also for {fj} = {0} the chain will be inhomogenous. 

4-3. Uniformization of the classical BS intertwining mapping. 

We like to establish an isospectrality transformation by commuting an auxiliary Lax 
operator through the monodromy Ai± f)109j) . For this we shall first parameterize the 
intertwining of two Lax operators considered in Sec. 13.11 maps (f79~J) and (jBHJ), in terms 
of cross ratios and if -functions. We use the results (I106J) specialized to a single site, e.g. 
n = 0. In this subsection this index will be suppressed. 
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For a compact notation, we introduce the following functions 

' j A-B H({fj(A) o-j(B)}) 

W({fj}, A, B) = - U ({fj} , B, A) . (110) 
For convenience we write A as the argument of U although it does not depend on A'. 
Using (jllUj) . equations (jl(J6|) uniformize the mapping (J75J) as follows: 

Ui, Wi = u, wtth}, y, z) ■ u 2 , w 2 = u, w({fj ^(y)}, x, z) ■ 

U 3 ,W 3 = U,W({f J }, x, y) ; 

U[,W{ = U,W({f j a j (X)},y, Z) ■ U' 2 ,W 2 = U, W({fj}, X, Z) ; 

u' 31 wi = u, w({fj o-jiz)}, x, y) , (in) 

with 

To do the same for the map (|83|) which describes the intertwining of the C, we 
go back to (|8D |) .(|8I ]) .(|52 |) and express these in terms of (JTTTJ). We shall write only the 
equations for the [/-variables, the W are analogous, see (|ll()j) . 

We describe the mapping TZ^ 23 (|HU|) by (jlllj) . (|112j) with the parameters 

{fj},x, y,z -. 

Ui = uttfj}, y, Z ) ; U 2 = U({f 3 a 3 (y)}, X, Z ) ; Us = Uttfj}, x, y) ; 

U* = U({f,a,(X)},y,Z ); U* = U({fj}, X , Zq) ; U 3 = U({fj <Tj(Z )}, X, y)- 

'—{S*}' K *-{H\' <-) 

The corresponding parameters for T^~ 3 will be called {gj} , X, y, and Z\\ 
U 1 = U{{g j },y,Z 1 ); U 2 = U({g j( r J (y)},X,Z 1 ) ; U^ = U ({ 9j } , X, y) ; 

Ul = U({g j <r j (X)},y,Z 1 ) ] U* = U {{ 9j } , X, Z,); U* = U ({ 9j a^)}, X ,y) ; 

R ^-{ Y z> l\ ^-{z[ l\ = (114) 

Since U'% of the first map (|113j) is the initial variable in the second map (|114|) we 
must take 

9j = fj o-j(Zo) . (115) 
Now the periodic condition (J52|) U 3 = U 3 requires that 

gj(Tj(Zx) = fj or o-j(Zo) o-j(Zi) = 1, (116) 
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which is equation (|92j) . Inserting (|115)) into (J 114)) we get the uniformization of the 
composite map s[% (JHSj) which solves the classical intertwining relations (fTSj) 

£ 2 (A)£i(A) = £*(A)£*(A). 
4-4- Auxiliary classical L-operator, Isospectrality 

We now introduce an auxiliary classical operator Cq ux (A) which by successive 
intertwining through the monodromy M. x ()1U9|) and imposing a periodic condition will 
lead to an isospectral transformation of the classical transfer matrix tr A^i(A). Using 
the notation (jllOj) the arguments of its operators Ci }n are 

Ui, n = u({fi utt °j\yi) },y n ,Zo), (H7) 

analogously Wi jn , see (jllOj) . U\,m W^i.n are obtained replacing Z$ — > Z\ and 

/, - /; ./>,(^o). 

We start writing (omitting the argument A which will always remain the same) 

paux p p r f* paux p p ( -i i 

^0 *-l,0 ••• *"1,Q-1 — *-l,0 ^1 • • ■ *-l,Q-l ■ 

We shall achieve our goal using as auxiliary operator the operator £2,7%, which is Ci )Tl 
of ()107|) with just the first indices "1" replaced by "2" (the second index ne 2 is not 
modified). Then the intertwining in ()118j) is the same as that considered in the previous 
subsection and the arguments of the various Lax operators are as follows (again, as in 
()117|h we write here only the Z7y ): 

£1,0 : U 1>0 = U ( {/,} , y , Zo ) ; C* 1>0 : £/£„ = U ( {f 3 a 3 {X)} , y , Z ) ; 

CT X : U 2fi = U ( {fMy )} ,X,Z ); CT X : U* >0 = U ( {/,•} ,X,Z ). (119) 

Moving C aux (A) one step further through the monodromy matrix (j 109)1 we write 
pis r>aux p p r r* r* paux p p 

*^l,0 *^1,2 • • • *-l,Q-l — ^1,0 ^1,1 ^2 *-l,2 • • • • 

Again using ()114j) and ()115)) . the variables of the new operators are: 

■ = U {{foWfaM} > Z o); CT X : Ul e2 = U ({fj/a^)} , X, Z ) . 
Continuing this way, we finally arrive at 

C% ux (A)Mi(A) = M\{A) C a $ x {A) , (120) 
where , Q1 . 

•^1 = ( II J C lo ^,2 • • • • • • £l,Q-l (121) 

V n'=0 / 

and the operators in (jl20|) . ()121)) have the following arguments: 

Ct,n ■■ Ul ne2 = U ({£ nr=i ^O*)}. 3>», 2b) ; (122) 

ar ■ ui ne2 = u ({£ nr=o ^W}. z ) . (123) 

We demand a periodic boundary condition for the auxiliary operator L a ^ x : 

C% X (A) = £« UX {A). (124) 
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This is fulfilled if the constraint (J91|) is imposed, i.e. if 

Q-l 

1,(3-1) = = i - 

Since the C^ x depend on the points y n only via {/j/o"j(iVn-i)} an d have the two 
other arguments fixed to be X, Z resp. X, Z\, the parameters K 2 , K 2 in the C°^ x are 
constant along the chain and are the same as in (jll3j) and (jll4j) . 

Taking the trace of (|12(Jj) and using (|124|) . we see that the transfer matrices 

TrA^i(A) and Tr M*(A) are isospectral. (125) 

A4i is composed of the Lax operators with the arguments ()117j) . the arguments of M\ 
are given in (|122|) . The classical integrals of motion Tj given by the generating series 
Q 

TjA' = Tr Mi{A) = Tr M$(A) (126) 

3=0 

are invariants of the isospectral transformation. 

Summarizing: equation (J92|) results from the periodicity in the e3-direction, while 
()91j) guarantees the periodicity of the auxiliary operator in the e2-direction. Both 
together determine the parameters Pj and P- . Due to the substitution fj — > f* = 
fj(Tj(X) in the variables, compare ()117j) with ()122|) . there is a non-trivial isospectrality 
if some fj ^ . 

4-5. Isospectral transformation of the BS quantum chain 

Finally, we generalize the isospectrality ()125j) of the classical BS-model to the quantum 
BS-model defined by the monodromy M(£) and transfer matrix T(£) of ()88|) . For this 
we have to find explicit expressions for the operators Q , L aux and M* in 

Q L aux (0 • M(0 = M*(0 ■ L* m (£) Q . (127) 

Once this has been found, imposing periodicity L aux (0 = £ aux (0 and taking the trace 
over both spaces C 2 and C , we will have the intertwining relation 

KT(0 = T*(£)K; T = Tr c2 M; K = Tr c ^ Q. (128) 

We start considering the intertwining of two quantum Lax operators ()89jl : 

Sjnn l J m\£i: "m) U m , . . . , K m ) • L n (J^ : 1t n , M n , . . . , fcn) 

= L n (£; it* , m*, . . . , k n ) ■ L m (^; u^n, . . . , k m ) S mn .(129) 

The Lax operators are matrices both in C 2 and C^. Written in components, (|129|) 
takes the same form as (|25jl . We also remark that in ()129j) . instead of L n of (jRHj) we 
could also use L(A; q n , q' n ) of (}2~4"j) since according to (jo^j) the gauge transformations 
of JUj) cancel: P~ 1 P m = 1 . 

We want to find an explicit expression for S mn for the case that the variables 
u m , u n , etc. on the left of ()129|) are related to those on the right, u^, w*, etc. , by the 
functional mapping (J8*3|) . i.e. by (jllHj) - pif))l . where we take m <->• 2, n <-> 1. 
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If the functional mapping is chosen to be trivial, (|129|) reduces to the Bazhanov- 
Stroganov intertwining relation (J25J) with S mn equal to S as given in (j2T?j) . (J4"2~)) . (J4"3~)) . 
Then in the cross-ratio parametrization the relation depends on the four pairs 
X, y, Zq, Zi . In the CP-parametrization one has to use (|68|). It appears that the 
parametrization of the intertwining operator S in terms of the CP-functions x p , x q , x p > 
and x q ' becomes inconvenient in the dynamical case since (jUHj) and (JoTJj) seem to have 
no simple generalization to the case {fj} ^ {0}. 

However, formulas (jfi7|) for the Fermat points have a nice dynamical extension. 
This can be used if we construct the quantum intertwiner S from the 3D-operator 
IZ123 imposing periodical boundary conditions. 

Now we can generalize the derivation of section l3~2l to the quantum case. In (JHJ) we 

have seen that since uj n = 1 , the map TZ\23 acting in the space of rational functions 

( f) 

of Weyl operators \l/ decomposes into the functional mapping 7Z\ 23 and a similarity 
transformation by the matrix R423 • The composition of two such maps is 



n i23 n 123 o $ = 7^ ( R ._ 3 . n m ( Ri23 . ^ . R -y . R -i 



123 

g> n[g o ( (n[{l o R i§3 ) • r 123 • * • ■ (n[{l o r^ 3 



(130) 



where the conjugation matrices map 

R123 : Ui.WnUa, w 2 ,u 3 , w 3 i-» u*, w*, u 2 , w 2 , u 3 , w 3 (131) 

and 

^123 oR -123 : Ui, W 1; U 2 , W 2 ,U' 3 , W3 ^ U^, W*,U*, W2,Ug, W* . (132) 

Imposing the periodic boundary conditions 

u 3 = u 3 ; W 3 = w 3 , (133) 

which imply the classical ones (JB2J), we define the quantum analog of the classical map 

~ by 



'12 



tr 3 (n$ 3 oR B3 ) -R 123 . (134) 



An expression (jZ 1 ^ o Rj 23 J means that the Fermat point parameters which enter the 
matrix elements of the matrix Ri 23 i n the third quantum space should be taken after 

( f\ - 

the functional mapping TCi 23 of ()14j) has been applied. 

Formula ()134|) provides a generalization of the Bazhanov-Stroganov intertwining 
operator (|42jl. It performs a canonical map of normalized Weyl operators 

4 = S 13 .^.SI?; § = 8x3.^.8^; etc. (135) 

It* Mi U{ U\ 

The Fermat points Xi, x 2 , x 3 determining the matrix Ri 23 in (|134|) are calculated from 
(JHHJ) raised to the iV-th power and then using (|113|) . Analogously for RJ23 take the 
iV-th power of (|H7jl and use (j!14jl . pi5jl . The result is 



jf _ j X Y'\ Hjjfjj^m H{{^{Y)a 3 {Z )}) 
Z< Z f F({f,(Y)}} //({f^Ay,^)^)}) 
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The R are determined by the Xj rather than by the x^ . So we must take iV-th roots. 
The possible choices of phases in this step have been discussed in |25] . 

The cancellation of the if-functions in the product x^x,^ arises as follows: The 

( f) 

matrix 7Z\ 23 °R-i23 has its Fermat parameters given by the same formulas as R123, just 
with Z replaced by Z\ , and the vector fj replaced by fj = fjO~j(Z ). Then take 
into account the periodicity o-j(Z ) Oj{Zx) = 1, equation (J9"2~j) . 

Given g, {fj} and the four pairs X, y, Z , Z\ , from ()136|) we obtain the Fermat 
points which are then used to calculate the operator S i2 from (jiSj) with (foTH). Suppressing 
to indicate g, Z , Z\ explicity, we shall write 

S 12 = S({fj},X,y). (137) 

Finally, we consider the BS-quantum chain of length Q. The Lax operators forming 
the monodromy are taken to be of the same form as the operator L n in fll29j) . with the 
scalar variables as in ()117j) . So, indicating the parameters, we write, recalling ()95|) : 

w 

^n(C; u n, u n , -z^, K n , k n ) = L(^; {fj Ijin - 1)} , y n , Z , Z x ). (138) 
so that the monodromy is 

m(0 = m{fj}, y ,z ,z 1 ) ■ L^,{f J /a J (y 1 )},y 1 ,z ,z 1 ) • ... 

■■■ • mifjljiQ-l^^Q^Z^Z^). (139) 

The auxiliary operator L aux is chosen to be of the same form, like L m in (j!29J) . with 
the scalar variables taken as for Cq ux in ()119)) : 

^(0 = L{t-{fjaj{y Q )},X,Z Q ,Z x ). 

Commuting the auxiliary operator through the chain proceeds in analogy to the classical 
case of section 14 A\ only in each step there is also the conjugation by a matrix 
S( {fjl(n — 1)} , X, y n ) . So M*(£) is given substituting on the right hand side of 
()139|) fj — > fj o~j(X) , and we have 

q = s({f 3 },x,y ) ■ sdfj/a^y^xM • s({f 3 i(Q-i)},x,y Q ^). (mo) 

Now all information necessary for the proof of Proposition El has been collected: 
The matrix conjugation by K does not change the spectrum so that we have to consider 
just the functional transform. Indeed, as seen in (jl()8|) u^Ui is independent of the 
{fj} and MiX has the same { /^-dependence as ttiX^ 1 and this is written in the first 
equation of (jHEJ)- Also, since Z n appears only in the combination (w n /w n ) Z n , from 
()108|) we confirm the second line of (}9l)|) . The two conditions (}9"T|) and (j92j) fixing the 
Pj, Pj have already been encountered in ()124|) and ()116|) . 

The operator K explicitly given by ()128)) . (jl4()j) . (jl34)) performs the isospectral 
transformation of the BS quantum transfer matrix since according to (j!28|) the spectrum 
of the transfer matrices T(£; { fj}) and T(£; { fj o~j(X) }) is the same. 
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5. Conclusion 

This paper has been devoted to describe the Bazhanov-Stroganov quantum chain using 
the tools of the 3D integrable generalized Zamolodchikov-Baxter-Bazhanov model in the 
vertex formulation of [Hj. The BS-Lax operator is constructed from the Linear Problem 
(J3J) imposing periodicity after two layers. The BS quantum intertwiner S , which is a 
product of four chiral Potts Boltzmann weights, is obtained applying twice the matrix 

conjugation part R 12 3 of the 3D mapping operator 7^123. The corresponding functional 

( f) 

operator 1Z 123 is used to solve the intertwining of two classical BS Lax-operators, a 
relation which would be difficult to obtain without the insight from 3D. 

There are many possible parametrization of the 3D mapping operator 7^123. These 
give rise to different more and less convenient parameterizations of the BS-transfer 
matrix. The parametrization in terms of cross-ratios and rational O-functions H turns 
out to be the most useful and is adopted for the explicit description of the isospectral 
transformations of the classical and quantum BS-transfer matrices. Whether our results 
form a sufficient preparation for solving the problem of separation of variables for the 
BS-model, similarly to what has been tried for other models in the papers (ElEHj, is 
subject of further work in progress. 
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Appendix: Alternative definitions of the Bazhanov-Stroganov L -operator. 

In this appendix we give the relation of two alternative formulations of the BS model 
to our definitions ()24|) . 

Appendix A.l. Bazhanov-Stroganov model 

In the Bazhanov-Stroganov original paper (B| the L-operator has been defined by 



where the Chiral Potts variables a p , b p , c p , d p , etc. are used, from which we obtain the 
variables of fT7jl . fTSj) bv 

x p = a p /d p ] y p = b p /c p ; ^ P = dp/c p , analogously with indices p, q, q' . 
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For the operators we use our notation X, Z of (0), so that XZ = cuZX.InjH] the 
Zjv operators are written X and Z, which correspond to our Z _1 and X -1 , respectively. 
Furthermore, p = (N — l)/2, and p\ is a constant. Extracting some factors, we get 



L BS (p; q,q') 
where 



p\uia p b p d q d q i Z p 



1 + AAs^Z 



A2X ( WIv — 



3/g 



A1X- 1 (x, - 
AiA 2 o;x (/ x„/ + 



Ai 



,1 




A simple gauge transformation by P = diag(X^ 2 , A 2 ^ 2 ) gives us the relation 

L BS (p; q,q') = p x uxi p b p d q d q , Z p PL(AiA 2 ; q, q') P' 1 
with our L(A; g, g') of (J23D • 



(A.3) 



Appendix A. 2. T^(t q ) -model 

In equation (5.33) of jH] a transfer matrix is introduced as the starting object for 
a fusion procedure, , which essentially is the transfer matrix of the BS model (|2T|) . 
has played a major role in solving the Chiral Potts model ^Hl UHl I2D1 121] ■ Since in many 
papers the definitions of [IB] have been used, we give the relation to our (|2~7j) . 
Equation (3.44) of [18 j is (as in (|2*7|). we denote the length of the system by Q ): 



U) 



k,q 



£ IF- 

mo,...,mQ_i=0 J=l 



2j(a 7+1 -V 7 +fc) Vq,j,<rj-<r'j+k 



T, 



Cj) 



Vq,j,mj 

= if j — k < a j — cr'j < N — k for any J. 



(A.4) 



The are transfer matrices leading from the initial Z N spins a = {a , . . . , <7q_i} to 
the final spins a' = {cr' , . . . , Cq_i} • There is a fusion hierarchy j = 0, . . . , N of the 
Tfc which is exploited in the applications to the CP. Of direct interest here is the case 
j = 2 , N > 2 . The index k labels a redundancy, see (3.51) of [TH], we shall take k — 0. 
Then the second part of (jA.4j) tells us that the Tq 2 J are non-vanishing only if for 

all J we have a j — a j = or 1 . So at J we can express the spin dependence in terms 
of the unit matrix and the two standard matrices Xj and Zj. 

Equation (3.37) of [THJ defines 

a+7V-j 

n (1-^) 



Vq,j,a 



with 



tq Xq l)q 



=a+l 



from which here we need only %2,i/%2,o = —ujt q . Equation (3.38) of [TH] gives 
F pg {j,a,m) = p%y; a - m ®{t p ,u a t q ) a J- a - x for < a, m < j < N, 
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where <&(x, y)™ 1 ' 17, is a polynomial in x, y which is expressible in terms of Gauss 
polynomials. We need only (equations (3.48) of [T5|): 

F M (2,0,0) = 1; F„{2^1) = -u% F w (2,l,0) = ^; F pq (2, 1, 1) = -J^v . 

y P y P y P 

Let us write 

1 i = 



-u;t 3 z = 1. 



Then, omitting the first argument j = 2 of _F(j, ...,...) and the index k = of g : 
[ T fL' = E w^W-^-^ F M ( S o,mo)F p , g ( Sl ,m Q )a;^K- 2 )^ FM ( Sl)mi ) x 

' n 1 ^Wln WWII 

mo,™i, m2...=0,l 

x F p , g (s 2 , m x ) a,"^""*) F M (s 2 , m 2 ) F^(s 3 , m 2 ) x . . . 

We collect all terms containing the spins o~j, o~j into a 2 x 2-matrix Tj : 

(rj) mj _ umj = uj m ^~ m ^^F pfq (sj,mj^)F pq (sj,mj) , (A.5) 

so that for periodic boundary conditions the transfer matrix is 

h (2) ] CT/j/ =Tr r T 1 T 2 ...T O _ 1 . (A.6) 
We now write ()A.5|) in matrix form: 

1 co^/y p 
-uj- aj ujt q / y p , -wt q / (y p y p/ 



( r ^)m. J _ 1 ,m J _,.,-<tj,.,# /,, , //',, „ A ^ CT J' CT J 



+ ~t q -co^x p § ^ 

y p y p > \ co aj ujx p it q co Sj cox p x p ' 1 ' ' 



mj-i,mj 



1/VV| 1 + A w Zj; - AX H^" ] z -i 

A ^-X/j^y - ^Z,) ; Au;* p ay + ^Z, 




' mj_ 1 ,m J 

where we put A = —l/{uot q ) , compare this with AiA 2 in (jA.2j) . So (L T denotes the 
2x2 matrix transpose of L ) 



\TJ ) m,j—i,mj 



1 AWV p jT p', p) P-%\ P= f iA " n 1/2 ° \<A.7) 



Observe that Baxter's variable £ 9 corresponds to our spectral parameter, hence the 
polynomial dependence of the transfer matrix on t q . 
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